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It is well known that the low-energy sector of quantum spin liquids and other magnetically disordered systems 
is governed by short-ranged resonating valence bonds (RVB). Here we show that the standard minimal trunca¬ 
tion to the nearest-neighbor valence bond (NNVB) basis fails completely even for systems where it should work 
the most, according to received wisdom. This paradigm shift is demonstrated for the quantum spin-1 square- 
kagome, where strong geometric frustration, similar to the kagome, prevents magnetic ordering down to zero 
temperature. The shortest tunneling events bear the strongest longer-range singlet fluctuations, leading to ampli¬ 
tudes that do not drop exponentially with the length of the loop L, and to an unexpected loop-six valence bond 
crystal (VBC), which is otherwise very high in energy at the minimal truncation level. The low-energy effective 
description gives in addition a clear example of correlated loop processes that depend not only on the type of 
the loop but also on its lattice embedding, a direct manifestation of the long-range nature of the virtual singlets. 


Introduction - The search for quantum spin liquids has been 
an active topic in condensed matter physics for many years [1- 
6, 18-21]. Recently, this search has gained new impetus from 
the discovery [7-15] of a series of layered spin-| kagome an- 
tiferromagnets (AFMs), whose ideal nearest-neighbor (NN) 
Heisenberg limit stands out as the prime candidate for a Z 2 
topological spin liquid [16, 17]. 

One of the earliest fundamental insights from the theory 
side is that the low-energy sector of all magnetically dis¬ 
ordered AFMs should be governed by resonances between 
sufficiently short-ranged singlet or valence bond (VB) pair¬ 
ings [6, 18-21]. Casting this very physical picture into an ef¬ 
fective Hamiltonian has been challenging for many years [22- 
31]. While one aspect of the problem, the non-orthogonality 
of the VBs [23], was essentially resolved recently [28-31], 
the really serious problem turns out to be the truncation to the 
NNVB basis, a problem that goes back to Zeng and Elser [24] . 

This is the second study that focuses on the impact of vir¬ 
tual singlets outside the NNVB basis to the tunneling physics. 
Ref. [31] has dealt with the 2D kagome, where the del¬ 
icate competition between different resonances leads to a 
VBC [32-35] which is very fragile [29, 30, 36], and virtual 
singlets are essential to stabilize the Z2 spin liquid [31]. 

The aim of this work is to show that virtual singlets have 
a strong qualitative effect in all disordered AFMs, even when 
the minimal NNVB truncation is expected to work the most, 
according to received wisdom. To this end, we have chosen 
an extreme case of an AFM, the square-kagome [37-42] of 
Fig. 1, which features a huge tunneling energy separation at 
the minimal truncation level. This happens because this sys¬ 
tem manifests the shortest resonances possible, the squares 
[□4 in Fig. 1 (b)], with NNVB amplitudes that are five times 
stronger than the second-shortest, Toop-six’ events. This en¬ 
ergy separation leads to the very robust ‘pinwheeT VBC of 
Fig. 1 (c, left) [37, 40], well separated in energy from compet¬ 
ing RVB states. 

Yet, even for such a system, virtual singlets change the 
physics qualitatively because they have very different impact 


on processes with different loop length L\ The square loops 
hold the largest density of defect triangles (triangles without 
VBs [43]) and thus bear much stronger quantum fluctuations, 
which suppress the loop-four amplitude by a striking factor 
of ^10. This leads to ‘renormalized’ amplitudes that do not 
drop exponentially with L, and to an unexpected Toop-six’ 
VBC [Fig. 1 (c, right)] which is otherwise very far in the pa¬ 
rameter space in the absence of virtual singlets. The result that 
virtual singlets change the physics qualitatively even in such 
a system simply tells us that the minimal NNVB truncation is 
probably inadequate for all standard [44] disordered AFMs. 

The ‘dressed’ NNVB description of the square-kagome of¬ 
fers in addition a clear example of correlated loop processes, 
which depend not only on the type of the loop but also on the 
particular VB environment of the loop. This embedding de¬ 
pendence is a direct manifestation of the long-range nature of 
virtual singlets, and adds another twist in the way we think 
about the effective RVB dynamics of disordered magnets. 

Model - The square-kagome features two symmetry in¬ 
equivalent sites (A and B), two NN bonds (AA and AB), and 
two void loops (AA-squares and AB-octagons), see Fig. 1(a). 
Hence the most natural spin-^ Heisenberg model comprises 
two NN exchange constants, Jaa and Jab, 

^Heis = ^ Jij^i • Sj , (1) 

Oj) 

where {ij) denotes NN sites. We target the nature of the disor¬ 
dered phase labeled ‘Y^’ in [40], which appears in the highly 
frustrated region x = ^ ^ 1. This is precisely where the RVB 
picture is most relevant, because at x = 1 we deal with a sys¬ 
tem of comer-sharing equilateral triangles, and such systems 
hold the two key RVB ingredients [24-26, 43]: i) a manifold 
of energetically favorable short-ranged VB states, and ii) a 
mechanism that drives the resonances between them. The first 
derives from the fact that each triangle can be minimized ex¬ 
actly by pairing two out of the three spins into a singlet. This 
leads to an extensive number of NNVB states where a large 
portion of the triangles (|) host a singlet and are thus satisfied 
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FIG. 1. (Color online) (a) The square-kagome has two inequivalent sites (A and B) and two NN exchange couplings Jaa and Jab- (b) The 
dominant processes around AA-squares. Blue and red ovals denote the two NNVB states involved in each process, (c) The two types of 
competing VBCs, the ‘pinwheeT (left) and the ‘loop-six’ (right), (d) Virtual excursions outside the NNVB basis. Here, the virtual singlet at 
the top is excited by applying the Hamiltonian on a defect triangle (shaded green) of one of the two NNVB states (red). 


locally. The second ingredient stems from the quantum fluc¬ 
tuations around the remaining, void or ‘defect’ triangles [43]. 

The square-kagome has 22 topologically distinct domi¬ 
nant processes, 4 around AA-squares and 18 around AB- 
octagons [40] . The former are the most crucial and are shown 
in Fig. 1 (b): The loop-four or ‘perfect square’ resonance 04 , 
two loop-six processes Dea and 066 , and Anally the loop-eight 
or ‘square pin wheel’ process Dg, involving NNVB states with 
2, 1 and 0 defect triangles, respectively. Similarly, the octago¬ 
nal processes [45] consist of: the ‘perfect octagon’ resonance 
Os. four loop -10 Oioa-d. eight loop -12 Oi 2 a-h^ four loop- 
14 Oi 4 a-d. and Anally the ‘octagonal pin wheel’ process Oie. 
with 4, 3, 2, 1 and 0 defect triangles, respectively. 

To accommodate the effect of virtual singlets on the quan¬ 
tum dimer model (QDM) [22-29, 50, 51] we write down an 
embedded QDM [31] 

£,e 

where we sum over all loops i that involve a resonance be¬ 
tween two NNVB states |1|) and |2|) [blue and red ovals in 
Fig. 1 (b)], and all possible VB environments e for given i. 
The corresponding tunneling amplitudes and potential ener¬ 
gies are denoted by and Vf (measured from Eo = — ^N,N 
is the number of sites). 

Minimal truncation & main competing phases - Within the 
minimal NNVB truncation method, the parameters and 
are independent of the embedding e [31]. Their values are 
reported in [40] based on a method which is equivalent to 
the inflnite order overlap expansion of [28, 29]. Namely, for 
each process we consider a minimal cluster that can host only 
the two NNVB coverings, |1|) and |2|), involved in the pro¬ 
cess, see top line of Fig. 2. For these minimal clusters, the 
NNVB problem amounts to diagonalizing the Hamiltonian 

^NNVB = ( . where O is the overlap ma- 

V / NNVB 

trix in the NNVB basis {11|), 12|)}. The matrix elements of O 
and Hnnvb can be found following standard rules [22-26,28]. 
The parameters involving the AA-squares are: = — 1, 

VOi = h ^Dea =*D6i. = = ^{2x-l), 

and tog = ^(1 — x), V\j^ — — x). For the remaining 

ones see Table I of [40]. 


Although the RVB description is only relevant around x = 
1 , it is useful to examine the phase diagram as a function of x, 
in order to establish the main competing phases of the prob¬ 
lem. Figure 3 shows exact diagonalization (ED) results from 
the 48-site cluster on a torus, using the parameters from the 
above 2x2 NNVB truncation. The low-energy spectra (a-b) 
and ground state (GS) expectation values of various loop den¬ 
sities (c) reveal three VBC phases: i) the ‘even pinwheeT of 
Fig. 1 (c, left) for X < 1, ii) the ‘odd pinwheeT of Fig. 1 (c, 
left) for 1 < X < 1.75, and iii) the ‘loop-six’ VBC of Fig. 1 (c, 
right) for X > 1.75. The potential terms do not alter this pic¬ 
ture qualitatively, as shown in Fig. 3. Without potential terms 
the transition between the ‘odd pinwheeT and the ‘loop-six’ 
crystals shifts to x ^ 2.34. 

All three states belong to the trivial topological sector, and 
break the translation symmetry in the same way, with a dou¬ 
bling of the unit cell. This symmetry breaking manifests in 
the spectra by an almost exact degeneracy between two states 
with momenta 0 and (tt, tt). The extremely small splitting [not 
visible in Fig. 3 (a,b)] shows that flnite-size effects are already 
very weak for the 48-site cluster. 

Fet us look at these states in more detail. The two ‘pin¬ 
wheeT states maximize the number of 04 resonances with 
even parity (i.e. with the symmetric combination of the two 
NNVB states involved, since = — 1 < 0 [52]), denoted 
by |PR+) in Fig. 1 (c, left). This amounts to a percentage 
of I resonating AA-squares. The remaining half carry ‘pin- 
wheels’, denoted by |PW±), whose parity is flxed by the sign 
of = ^(1 — x), i.e., it is odd (even) for x < 1 (> 1 ). 
At X = 1 , the two parities give rise to an extensive degener¬ 
acy of 212 (on top of the two-fold degeneracy from transla¬ 
tional symmetry breaking), which can be clearly seen in the 
exact spectra of Fig. 3 (a,b). This is analogous to the 36-site 
VBC of the kagome [30, 32-35], which maximizes the ‘per¬ 
fect hexagon’ resonances, and features a similar extensive de¬ 
generacy due to ‘hexagonal pinwheels’. Here the degeneracy 
is lifted for X 7 ^ 1 due to the inequivalent bond structure [40]. 

Turning to the ‘loop-six’ crystal, this state maximizes the 
□ 6 a resonances, denoted by |L 6 _) in Fig. 1 (c, right). Here all 
AA-squares carry such a resonance and the parity is odd since 
taea = I ^ 0. We emphasize here that at the level of the min¬ 
imal NNVB truncation, foe favored ‘loop- 
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FIG. 2. (Color online) Finite spin-1/2 Heisenberg clusters from which we extract the tunneling amplitudes for the four dominant processes on 
the lattice, see text. The shaded (green) triangles indicate the closest extra defect triangles that give rise to an embedding dependence. 


six’ phase features only resonances of the type Dea- The latter 
are therefore selected by weak quantum fluctuations driven by 
the remaining tunneling terms in the minimal QDM. Indeed, 
states featuring Deb resonances or combinations of Dea and 
□66 have a very small excitation energy, as can be seen in the 
low-energy excitation spectrum of Fig. 3 (a,b). 

Fluctuations are in fact very weak for all three VBC phases. 
As shown in Fig. 3 (c), the GS loop densities of 04 and Ds are 
almost exactly equal to 1/2 for x < 1.75, while that of Dea is 
almost exactly equal to 1 for x > 1.75. The same fact is also 





FIG. 3. (Color online) ED results of (2) using the minimal NNVB 
parameters, as a function of x = Jab / Jaa, with and without potential 
terms V'. (a-b) Low-energy spectra measured from the GS. The 32 

N 

states (2 X 212 ) that become degenerate at x = 1 (see text) are shown 
in grey, and the number d is the degeneracy of each level, (c) GS ex¬ 
pectation values of various loop densities, (d) Comparison between 
the GS energy with (3), see text. 


demonstrated in Fig. 3 (d) which shows that the exact energies 
of the lowest eigenstates (symbols) lie practically on top of the 
corresponding energies (lines) 

<°nlheel/^ = ” ^4 I + Vq, - ItaJ) /12 , 

4op-six/^=(^a--|^a-l)/6’ (3) 

of the ideal product states between the resonating building 
blocks without extra quantum fluctuations. The latter are very 
weak because they arise from sequences of longer loop events 
that necessarily involve the octagons, whose tunneling ampli¬ 
tudes are very small, see [40] . 

'Dressed’ NNVB description - Having established the main 
competing phases of the problem we now focus again on the 
highly frustrated point x = 1. As announced above, the longer- 
range singlets that are virtually excited around the defect tri¬ 
angles destabilize the ‘pinwheel’ crystals in favor of the res¬ 
onating Toop-six’ state, i.e. they effectively shift the boundary 
between the two phases in a drastic way. Essentially, the ratio 
~ta 4 /ta 6 a drops from the NNVB value of 5 below 2, which 
is the transition value obtained from (3) with and by 

neglecting the potential terms (see below). 

Let us see how such a drastic reduction takes place here. 
Following [31], we extract the QDM parameters from the ED 
spectra of appropriate flnite Heisenberg clusters, using the 
lowest two exact eigenvalues with the right symmetry. This 
method is essentially equivalent to enlarging the truncation 
basis in order to include the virtual singlets accommodated by 
the clusters. Figure 2 shows the clusters considered here for 
the four processes around AA-squares (for the octagon pro¬ 
cesses see [45]). The clusters shown at the top line capture 
the virtual singlets that are closest to the central loops. 

To capture the contribution from other virtual singlets we 
proceed as in [31] and add sawtooth chains of length R at the 
positions of the defect triangles of the minimal R = 0 clusters. 
These chains provide the only escaping paths of the virtual 
singlets (or the spinons forming these singlets [53]), provided 
we do not encounter extra defect triangles nearby. Such ex¬ 
tra defect triangles become first relevant for the loop-six pro- 
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R=0 

R=1 

R=2 

R=3 

□ 4 

-1 

-1 

-0.385403 

-0.208442 

-0.133124 

^6a 

+1/5 

+0.266044 

+0.151461 

+0.125542 

+0.113368 

^6a 

+1/5 

+0.266044 

+0.131499 

+0.102100 

+0.089369 

^6a 

+1/5 

+0.266044 

+0.124130 

+0.093471 

+0.081068 


+1/5 

+0.116951 

+0.097390 

+0.088453 

+0.082762 

^6b 

+1/5 

+0.116951 

+0.063010 

+0.057934 

+0.054884 

^6b 

+1/5 

+0.116951 

+0.002895 

+0.029735 

+0.030038 

□ 8 

0 

0 

0 

0 

0 


TABLE L Tunneling parameters t for the square processes (for the 
octagonal processes see [45]), as extracted from the finite spin-1/2 
Heisenberg clusters of Fig. 2, as described in the text. The numbers 
from the minimal NNVB truncation are also shown for comparison. 
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FIG. 4. (Color online) Connected dimer-dimer correlation profiles in 
the GS of (2) at X = 1, using the NNVB tunneling parameters (left) 
and the ones extracted from the i? = 2 clusters of Fig. 2 (right), which 
includes the explicit dependence on the extra defect triangles (see 
text). The reference dimer is the horizontal segment at the bottom 
left corner. The red and grey (left) or blue and grey (right) denote 
positive and negative values. 


cesses, and can be studied by examining different variations 
of the corresponding clusters for each given R, that accom¬ 
modate zero, one or two extra defect triangles, see Fig. 2. The 
corresponding effect on the octagonal processes is disregarded 
here since these processes are anyway too weak to affect the 
low-energy physics quantitatively. 

The extracted tunneling amplitudes for the square and the 
octagonal processes are given in Table I and Table 2 of [45], 
respectively, and deliver a number of clear insights (the po¬ 
tential terms are much smaller and can be disregarded [45]). 
First, except for the minimal = 0 embedding of (and 
the processes Dg and Oie whose amplitudes are zero irrespec¬ 
tive of R), the amplitudes are clearly different from the ones 
resulting from the minimal 2x2 truncation (NNVB). The dif¬ 
ferences show that the 2x2 truncation give in most cases a very 
poor approximation, so the virtual singlets cannot be ignored. 

Second, the shortest amplitude shows a remarkable re¬ 
duction by about ten times from R = 0 to R = 3. By contrast, 
the amplitudes drop by about 2-3 times (depending on e), 
and show a relatively faster convergence with R compared to 
Virtual singlets have therefore a much stronger tendency 
to penalize the shortest loop events because they correspond 
to the highest local density of defect triangles. With the ra¬ 
tio dropping below the transition value of 2 above 

R = 2, it is immediately clear that the RVB physics will no 
longer be governed by the shortest resonances. 

Third, our results show a ~ 20% variation of and a 
much stronger ^ 40-60% variation of depending on e, 
i.e. on whether there are zero, one or two extra defect tri¬ 
angles nearby. Such a strong dependence on one class of the 
dominant processes of the problem must be taken into account 
explicitly in the effective description. The presence of such 
correlated loop processes, which depend strongly not only on 
the type of the loop but also on the embedding, highlights the 
qualitatively non-trivial impact of virtual singlets, which ef¬ 
fectively sample the immediate VB environment of a loop. 

Fourth, virtual singlets have a different impact on and 
again in contrast to the minimal NNVB truncation. For 

(e) (e) 

ah R we systematically find that showing a clear 


preference for loop-six resonances of the u^a typ^- In the 
minimal NNVB model discussed above, such events were se¬ 
lected by weak virtual octagonal processes whereas here they 
are favored directly by virtual singlets. 

Finally, the processes around octagons are also renormal¬ 
ized substantially by virtual singlets, see [45]. But, these pro¬ 
cesses remain anyway harmless, which is checked by a full 
ED of the ‘dressed’ QDM including ah 22 processes. 

Fig. 4 shows the connected dimer-dimer (or density-density 
in the QDM setting) GS correlation profiles using the param¬ 
eters from the minimal NNVB truncation and from the R = 2 
clusters, with the explicit embedding dependence for the loop- 
six processes. At this level, we have access to ah amplitudes 
except /q , whose corresponding cluster is too large for ED. 
In the calculations shown here we have taken the conservative 
assumption which does not affect the results. 

The data confirm that the ‘loop-six’ crystal is stabilized al¬ 
ready at R = 2, and taking the i? = 3 amplitudes renders this 
phase even more robust. 

Outlook - We have shown that the minimal truncation to 
the NNVB basis fails qualitatively even in a system where it 
is expected to work the most, according to received wisdom. 
This is strong evidence that virtual singlets have a qualitative 
impact in ah standard disordered AFMs. Not only do the vir¬ 
tual singlets ‘dress’ the basic RVB building blocks, they also 
completely alter fundamental notions that derive from our ex¬ 
perience with the minimal NNVB truncation: The tunneling 
amplitudes do not drop exponentially in the length L of the 
loops, and the tunneling dynamics cannot always be approxi¬ 
mated as a sum over strictly local loop processes, as initially 
envisaged by Rokhsar and Kivelson [23]. 
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Supplemental material 

A. General remarks for the NNVB manifold of the square-kagome antiferromagnet 

As shown in [1] by a generalization of the so-called arrow representation [2], the dimension of the NNVB Hilbert space is 
2 X 2^/^, where N is the total number of sites. Similar to the non-bipartite, triangular and kagome lattices, the NNVB manifold 
splits into four topological sectors for periodic boundary conditions in both directions. These sectors can be labeled by the 
parities (Px,Py) of the number of dimers intersecting any horizontal or vertical line [3-6]. Here, the sectors (0,1) and (1,0) are 
connected to each other by the four-fold symmetry of the lattice. 


B. Basic elements of the minimal 2x2 NNVB truncation 


The basic elements of the minimal 2x2 NNVB truncation method in the square-kagome AFM are given in [1]. For convenience 
they are repeated here in Table II. Let us consider any of the finite spin-1/2 Heisenberg clusters of the first row of Fig. (2) of 
the main text (or Fig. 5 for the octagonal processes) and consider the elementary tunneling process p between the two dimer 
coverings \lp) and \2p) that can be accommodated by this cluster. For the dimer orientations [7], we take the singlets to be 
oriented clockwise in each void square and octagon. The transition graph [8] of our elementary process involves a single non¬ 
trivial loop of length Lp, surrounding a single square or a single octagon. Then, within our convention, the overlap between the 
two NNVB states is given by 


= ( 1 ,| 2 ,) = (- 1 )^ 21 -^ 


The diagonal and off-diagonal matrix elements of the Hamiltonian, 

Eo^p = (lp|7/Heis|lp) = (2p|7/Heis|2p), Vp = (lp|7/Heis|2p), 
can be found using standard rules [8-14]. In the basis {|lp), |2p)} we have 


with 


^NNVB = (O ^ )nnvB = 


tp ~\~hr 


1 Ldp 

LUnj 1 


^0,p 

^p Eo,p 


1 Ldp 

Ldnj 1 


^0,p 3 “ ^p tp 

tp ^0,p “1“ ^p 


l-ul 


, Vp = —tp udp, where hp = Vpjudp — Eq^p . 


( 4 ) 

( 5 ) 

( 6 ) 

( 7 ) 


The quantity hp corresponds to Eq. (12) of Schwandt et al[\5], modulo a factor of 3/4 related to the redefinition of the Hamilto¬ 
nian; see passage below Eq. (10) in [15]. So Eq. (7) is consistent with Eq. (40) of [15]. In other words, the above 2x2 truncation 
is equivalent with the infinite-order cluster expansion of [15]. 


void plaquette 
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TABLE II. QDM parameters tp and Vp (in units of Jaa) for the four processes around AA-squares and the 18 processes around AB-octagons, 
using the minimal 2x2 NNVB truncation approach. Lp is the length of the loop in the transition graph of |lp) and \2p). For the dimer 
orientations we follow the convention that singlets are oriented clockwise in each even-length loop. All values correspond to the clusters of the 
first row of Fig. 5. If we use clusters of any other row, Eo^p and Vp will change, but the overlap cOp and the combination hp = Vpjujp — 
do not change, leading to the same tp and Vp. 
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C. Potential terms 

Here we discuss the potential terms of the effective QDM. As explained in [16], the largest contributions to the potential 
energy arise from processes involving a single defect triangle, and as such they can be absorbed in a global energy shift, since 
the total number of defect triangles is a constant in the NNVB basis. The remaining contributions can be divided into binding 
energies among two, three, etc. defect triangles. In the present case, the most important processes (i.e. the ones around the AA- 
squares) involve up to two nearby defect triangles, for which the binding energy is about 0.01 [16], i.e. 10 times smaller than the 
dominant tunneling amplitudes. This is why we can safely disregard the potential terms in the effective QDM description of the 
square-kagome. 


D. Finite clusters for the octagonal processes 

Figure 5 shows the finite-size spin-| Heisenberg clusters that are used to extract the tunneling parameters of the 18 different 
octagonal processes (for the processes around AA-squares see main text). 
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FIG. 5. Finite clusters used to extract the octagonal tunneling parameters of the effective Hamiltonian. The gray triangles are the defect 
triangles involved in each tunneling process. 








